Rise and fall of hidden string order of lattice bosons 
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We investigate the ground state properties of a newly discovered phase of one dimensional lattice 
bosons with extended interactions^. The new phase, termed the Haldane Insulator (HI) in analogy 
with the gapped phase of spin-1 chains, is characterized by a non local order parameter, which can 
only be written as an infinite string in terms of the bosonic densities. We show that the string 
order can nevertheless be probed with physical fields that couple locally via the effect those fields 
have on the quantum phase transitions separating the exotic phase from the conventional Mott 
and density wave phases. Using a field theoretical analysis we show that a perturbation which 
breaks lattice inversion symmetry gaps the critical point separating the Mott and Haldane phases 
and eliminates the sharp distinction between them. This is remarkable given that neither of these 
phases involves broken inversion symmetry. We also investigate the evolution of the phase diagram 
with the tunable coupling between parallel chains in an optical lattice setup. We find that inter- 
chain tunneling destroys the direct phase transition between the Mott and Haldane insulators by 
establishing an intermediate superfluid phase. On the other hand coupling the chains only by 
weak repulsive interactions does not modify the structure of the phase diagram. The theoretical 
predictions are confirmed with numerical calculations using the Density Matrix Renormalization 
Group (DMRG). 



I. INTRODUCTION 

Systems of ultracold bosons in optical lattices offer 
unique opportunities for studying strongly correlated 
quantum matter in a highly controllable environment 1 . 
In a sufficiently deep lattice potential the interactions 
dominate over the kinetic energy and, at commensurate 
filling, can drive a quantum phase transition from the 
superfluid to the Mott insulating state. This transition 
has been observed experimentally 2 -, and it requires only 
local (on-site) interactions between the bosons. Longer 
range interactions, effective for example in systems of ul- 
tracold polar molecules^ or atoms with large magnetic 
dipolc moments* , can give rise to even richer behav- 
ior. In particular, we have recently predicted that bosons 
with sufficiently strong nearest neighbor or further range 
repulsion on a one dimensional lattice form a new insu- 
lating ground state characterized by hidden topological 
order—. 

We termed the new phase a Haldane Insulator (HI) be- 
cause of the close analogy with Haldane's gapped phase 
of integer spin-chains^. Both states support a highly non 
local string order paramete r 8 ! 9 ' 10 . Moreover, the phase 
transitions from the Haldane insulating phase to the con- 
ventional Mott and Density wave (DW) insulators, also 
have their analogies in anisotropic spin chains 8 - ! 12 ' 13 . This 
phase can potentially be realized in a system of cold dipo- 
lar atoms or molecules in a one dimensional optical lat- 
tice, where the dipole moment is polarized perpendicular 
to the chain direction. 

There are, however, essential differences between the 
spin chains and the lattice bosons. First, the anisotropic 
spin chains enjoy a global Z2 symmetry associated with 
flipping the S z component of all spins. This would trans- 



late to a particle-hole symmetry about the mean lattice 
filling, which is clearly absent in the microscopic bosonic 
models. Furthermore, the optical lattice systems are in 
practice not strictly one dimensional. Some amount of 
residual coupling between chains due to tunneling and, 
in our case, also dipolar interactions, is inevitable. More 
generally, we shall see that additional external fields, 
which can be applied to the ultracold Bose systems, cou- 
ple in a non trivial way to the string order parameter. 
The effects of these perturbations on the phase diagram 
raise new fundamental questions on the nature of the non 
local order. 

In this paper we investigate how different perturba- 
tions, natural to bosonic systems, affect the transitions 
from the Haldane insulator to the conventional Mott and 
Density wave phases. We argue that these directly probe 
the nature of the non local order parameter which charac- 
terizes the Haldane insulator phase. Using a bosonization 
analysis, we demonstrate that the transitions between the 
HI and the conventional phases are not sensitive to the 
broken particle-hole symmetry in the microscopic Hamil- 
tonian. Such perturbations also leave the string order 
parameter intact. On the other hand we predict that per- 
turbations that break the lattice inversion symmetry in 
addition to the particle-hole symmetry, would eliminate 
the distinction between the MI and HI phases and gap 
out the critical point separating them. Next we inves- 
tigate the effect of coupling between two parallel chains 
("ladder" geometry). The predictions of the field theoret- 
ical analysis are confirmed with numerical simulations us- 
ing the Density Matrix Renormalization Group (DMRG) 
method. 

Our starting point for theoretical analysis is the ex- 
tended Bose Hubbard model (EBHM) on one or two 
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FIG. 1: (Color online.) Array of one dimensional chains 
formed by a two dimensional optical lattice. The polariz- 
ing field E is directed perpendicular to the chains and at an 
angle which can be used to tune the inter-chain interaction 
V ± in ©. 

chains at single site occupation (n = 1): 

3 

+^1]^ ( Ha i _ X ) + V J2 n ocjn aj+l , (1) 

3 3 

where j is the site index and a = 1 , 2 is the chain index, 
in the case of two chains. The inter-chain coupling in 
that case includes a transverse tunneling matrix element 
and interaction term: 

Hx = -t x $>i>i + H.c.) +VxJ2 

n lj n 2 j. (2) 

3 3 

In practice the intcr-chain tunneling t± may be tuned 
by varying the depth of the optical lattice potential in 
the direction perpendicular to the chains. The inter- 
chain dipolar interaction can be controlled by changing 
the direction of the external polarizing field relative to 
the plane of the chains, while keeping it perpendicular to 
the chain axis. A schematic setup of this type involving 
many coupled chains is illustrated in Fig. [T] The two 
chain model 

H= ]T H a + H ± (3) 

a=l,2 

serves as a bridge between one and two dimensional ge- 
ometries, while still amenable to powerful techniques for 
treating one dimensional systems, such as bosonization 
and the density matrix renormalization group (DMRG). 
Similar quasi one-dimensional systems were studied pre- 
viously in different parameter regimes: with only on- 
site interaction s 14 ' 15 and with longer-range interactions 
at half-integer filling^. However, the Haldanc insulator 
phase which is the focus of this paper is not obtained in 
these regimes. 

Recently Anfuso and Rosch noted that the string or- 
der characterizing the Haldane gapped phase of spin-1 
chains is unstable against inter-chain antiferromagnetic 
exchange^. The string order in the HI phase is similarly 
sensitive to inter-chain tunnel coupling (see appendix [A"]) , 
and the distinction between the HI and MI phases is lost. 
What is then the fate of the direct second order transi- 
tion found between these two phases when intcr-chain 



tunneling is turned on? Interestingly we find that di- 
rect transition between these two phases is avoided in 
the double chain system by the appearance of an inter- 
mediate supcrfluid phase for arbitrarily small inter-chain 
tunneling. 

Our analysis consists of the following parts. In section 
III Al we present a low energy effective field theory for the 
single chain, which correctly captures all three insulating 
phases MI, HI and DW. This is done with help of a spe- 
cial bosonization procedure borrowed from work on inte- 
ger spin chains 1 ^. The bosonized forms for the non local 
correlations that characterize the HI and MI phases are 
given in section Til Bl This framework is used in section 
III CI to study the coupling of the non local order parame- 
ters to inversion symmetry breaking perturbations. Then 
in section III Dl we extend the bosonization approach to 
describe a pair of weakly coupled chains. We carry out a 
renormalization group (RG) analysis to obtain analytic 
predictions for the phases arising at weak inter-chain cou- 
pling in the vicinity of the quantum critical points. The 
most interesting result of the coupling is the effect of 
inter-chain tunneling on the transition between the MI 
and HI states. This perturbation is highly relevant at 
the critical point. As a result we find that for any fi- 
nite tx the HI and MI insulating phases are separated 
by a superfluid region, whose domain grows rapidly with 
increasing t± . Finally in section IIIII we confirm the ana- 
lytic predictions and extend them to stronger inter-chain 
coupling using numerical DMRG calculations. 

The paper is followed by three appendices. In appendix 
[A] we investigate the effect of weak inter-chain coupling 
deep in the HI phase within perturbation theory. Al- 
though the thermodynamic HI phase is protected by a 
gap A, we show that the single chain string order is de- 
stroyed even by infinitesimal coupling. In appendix [B] 
we provide an alternative derivation of the effective field 
theory using direct bosonization of the particles rather 
than resorting to an effective spin-1 model as done in 
section [H] In appendix [U] wc construct explicit expres- 
sions for the string order parameter within the effective 
field theory. 



II. BOSONIZATION 
A. Continuum limit of a single unperturbed chain 

A single chain of interacting bosons |T]) at filling ?i = 1 
was studied in Ref. Q using DMRG and the existence 
of the new Haldane insulator (HI) phase was predicted. 
The phase diagram in the two dimensional space of U/t 
versus V/t is reproduced here in Fig. H{a). 

The insulating phases and the correct quantum phase 
transitions separating them can also be obtained from a 
field theoretical analysis^, which builds on a direct anal- 
ogy with the bosonization procedure developed for inte- 
ger spin chain s 11 ! 12 ' 20 . Let us briefly review the deriva- 
tion of the effective field theory of a single chain. The 
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first step is an approximate mapping of the EBHM ([T]) 
to an anisotropic spin-1 model. This is done by trunca- 
tion of the Hilbcrt space of each site to the three lowest 
occupation states n = 0, 1, 2 21 i 22 . Such a truncation is 
justified at large U when fluctuations in the site occupa- 
tions are strongly suppressed. The effective spin-1 model 



Sj~S j+1 



IE 

3 



3+1 



(4) 

Here S z = n — n. In general, there are other terms in 
this Hamiltonian which break the S z — > — S z , S + < — > 
S~ ( "particle- hole" ) symmetry. However, in the effective 
long-wavelength Hamiltonian we will now derive, these 
terms are irrelevant, so we will drop them for now. The 
effect of these terms will be considered in section III CI 
One then proceeds by writing each spin 1 variable as 

Each 



1,3 



'2,3' 



a sum of two spin | variable o n i 12 , S" = s 
spin-i chain can now be mapped to a spinless fermion 
chain by a Jordan- Wigncr transformation. The fermions 
are then bosonized in the standard w&yH, according to 
i>a,R/L ~ v /27ra e ^ e ° ± ^°' > ' wnere a * s tne lattice spac- 
ing, a = 1,2 corresponds to the two fictitious spin-1/2 
chains (not to be confused with the two physical chains 
that will be considered in section HlD|) and R/L corre- 
sponds to right and left moving fermions, respectively. 
The bosonic fields satisfy the canonical commutation re- 
lations [4> a (x) , Op (x 1 )] = —iirSa^Q (x 1 — x), where O (x) 
is a Heaviside step function. With these conventions 2 ^, 
the spin-1/2 operators have the following bosonized form: 



-d x 



(-1) 1 



sin (2(j) a ) 



(5) 



The bosonized Hamiltonian is conveniently written in 
terms of symmetric (+) and antisymmetric (— ) combi- 
nations of the fields <f) a and a : 



H = H+ + H-+H+-. 
The three terms are given explicitly by: 



(0) 



Ztt 



K + (d x + f + -L(d x + ) 2 



dx 



9i 



H- = 



2tt 



(no) 
dx 



■ cos (2</) + ) 



(7) 



dx 



^2 cos (2<£_) + j^-s cos {20 J) 1(8) 



(tto) 



and 



where 



H 4 



dx- 
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(tto) 
hi ± (f>2 and 0± 



cos (2(j)+) cos {2(f)-) 



(9) 



h ± ^2) /2. (Note that 



the fields </> + and 0_ correspond to v / 2V'i an d V2ip2 in 



Ref. [13, respectively.) The naive continuum limit gives 
the following estimates for bare values of the coupling 
constants: 



U + 6V T ^ 
u + = ta\jl-\ — — , A 4 



U+6V 

7Tt 



U-2V T , 
u_ = ta\\\ , K _ = 

7rf 



U-2V 

TTt 



(2V-U)a 
9i = -.92 = 7, 1 33 



—tira, 34 = Va(10) 



To study the physical correlation functions we need 
to express the second quantized Bose operators on the 
lattice and the local site occupation in terms of the con- 
tinuum fields. These relations arc of the form 



Hz) 

n (x) 



-d x 4>+ 

7T 



COS ( 



-) + 



(11) 



-) cos (</>_) + 



where A and B are non-universal constants and only the 
most relevant terms are shown. In general, the expan- 
sion of these operators also contains less relevant (sub- 
leading) terms. 

The Hamiltonian H+ of the symmetric degrees of free- 
dom is similar to the usual low energy description of lat- 
tice bosons. The Umklapp term parameterized by gi 
is relevant for K + < 2. Here we consider the insulating 
phases, for which K + is indeed below this critical valued. 
The main difference from the standard bosonization of 
lattice bosons is that the parameter g\ changes sign at 
2V = U. This change of sign marks a quantum phase 
transition from the MI, for which + is localized around 
zero, to another gapped phase in which it is localized 
near ±~. As we shall see below, the latter turns out to 
be the Hi phase and is characterized by string order. On 
the critical line, the system is described by a Luttingcr 
liquid with power law correlations that depend on the 
non-universal parameter | < K + < 2. If K + is driven 
below i, the critical line becomes unstable, since the sub- 
leading term cos(4</> + ) (which was omitted in Eq. 0) 
becomes relevant. This possibility will be discussed fur- 
ther in Sec. Ill CI Note that the odd sector described by 
is not critical. The term g% in © is relevant on both 
sides of the transition, so that the field 0- is localized 
and remains massive at the critical point. 

On the other hand, the transition from the HI to the 
DW phase is controlled by the odd sector. It occurs when 
the term cos (2(f)-) becomes more relevant than the dual 
term cos(26L), that is when K- < 1. According to 
the naive continuum limit the critical line is given by 
V — 2U = 3irt. This transition is of the Ising universality 
clas o 12 i 24 . This can also bee seen by refermionization of 
the Hamiltonian ii_ at the critical point K_ = 1 and 
.92 = .93, which yields precisely the masslcss Majorana 
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fermions of the 1 + 1 dimensional Ising model at critical- 
ly. 

We remark that the precise values of the coupling con- 
stants as a function of the microscopic parameters are not 
simple to derive. The field theory ^ should be used to 
obtain the universal behavior of the system near the crit- 
ical points which describe the quantum phase transitions, 
not to find the precise location of phase boundaries. 



B. The string and parity order parameters 

What is the physical distinction between the different 
insulating phases? It was shown in Ref. Q that the HI 
phase of bosons is characterized by a non decaying string 
correlation function 



0\ 



lim ( 5rii exp | iir y Sn^ | Srij 



i<k<j 



(12) 



Here Srii = rii — n is the deviation from the average (in- 
teger) filling n. In section Hill we show that the MI phase 
is characterized by a different non local correlation func- 
tion, which will be termed the "parity" correlation func- 
tion 



Oi> = lim ( exp iir > 8ni 



(13) 



For the sake of convenience, we define the parity and 
string operators 



A simple interpretation of these correlations is illus- 
trated in Fig. [2J In both phases the ground state con- 
sists of configurations with most sites having precisely 
the average occupation, but also some particle and hole 
fluctuations. The parity Op order in the MI phase im- 
plies that the fluctuations in this ground state are bound 
particle-hole pairs. On the other hand, the string order in 
the HI phase implies that a "renormalizcd" chain with all 
non fluctuating 5rii — sites taken away would have den- 
sity wave order (particle and hole fluctuations alternating 
along the chain). Physically, this can be understood as a 
compromises between the kinetic energy term in Eq. ([1} , 
which prefers maximum derealization of particles and 
holes, and the nearest neighbor interaction term, which 
is minimized when particles and holes are neighbors. In- 
deed, the HI phase is realized at intermediate values of 
V/t. 

Within the effective field theory the HI and MI seem 
to differ only in the expectation value of the field </>+ 
which is pinned in each of these phases. It is interesting 
to relate this distinction to the string and parity corre- 
lations that characterize the two phases. Since both or- 
der parameters contain the factor exp (ift^n^Kj ^ n kj , 
we may naively expect their bosonized forms to con 
tain exp [i 

be symmetrized carefully to obtain a hermitian opera- 
tor. In Appendix [Cj we argue for the following forms of 
Os and Op in the bosonized theory: 



+ 0)], since T,i<k<j 6n k ~> h fx* dxd *<t>+ = 
- 4>+ [xi)\. However, the exponential should 



O s {x) 



(<t>+ (*)) 



(16) 



&p u) = cx p ( in 

k<j 



(14) 



Op(x) ~ cos ((/>+ (x)) 



(17) 



6 s (j) = 6 P {j)8n j , (15) 

such that 2 P = lim^^ (dp (i)6 P and 0% = 
lim K _ iHoo i&s (i)O s (i)}. 

+ + 0)®® ® - I 

(a) 

+ - (0) ® ®0 ® (±)}© © © (0-- I 

(b) 

FIG. 2: Typical configurations in the HI (a) and MI (b) 
ground states. The numbers represent 8n (the deviation of 
the local occupation from the average density). The HI can 
be described as a charge ordered +,—,+,... state with an un- 
determined number of sites between each + and — . The MI 
is a dilute gas of particle-hole pairs (indicated by the dotted 
line) . 



Here the form of ([To]) was postulated on the basis of 
symmetry. For a more microscopic derivation see Rcfs. 

Sin. 

The above expressions for the string and parity cor- 
relations are consistent with the respective phases de- 
rived from Eq. (|7l8p . The MI phase corresponds to 
gi < 0, which implies a non vanishing expectation value 
of the parity operator (cos (<f) + {x))) ^ 0. The fact 
that (/>+ is locked to or 7r in this phase is consistent 
with the particle density being concentrated on the lat- 
tice sites, as in the cartoon product state | ^mi ) ~ 
fjj -bj | ) . In the HI phase on the other hand g\ > 0, so 
that(sin (</>+ (x))) ^ and therefore non vanishing string 
order. Here <f) + is locked to 7r/2 or 37r/2, which implies a 
shift of the particles by half a lattice constant compared 
with the MI, that is, the density is centered on the links 
rather than the lattice sites. This is captured by the car- 
toon wave function | ^ HI ) ~ rij(^] + I ® ) > which 
is closely analogous to the AKLT state that describes a 
valence bond solid in spin-1 chains. 
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C. Coupling to symmetry breaking perturbations 

A realization of the phases described above with ultra- 
cold atoms would open up new ways to probe the nature 
of string and parity orders by how they react to different 
perturbations. Interestingly, we shall see that in spite 
of being highly non local, the string and parity operators 
nevertheless couple in interesting ways to local symmetry 
breaking fields. 

As compared to the spin chain model [Eq. , an in- 
herent broken symmetry in the microscopic Hamiltonian 
([1]) is the absence of particle hole symmetry Srii — > —Srii. 
The spin model does enjoy the analogous symmetry as- 
sociated with rotation by tt around the y axis, which 
takes S z — > — S z and S x — > — S x . Nevertheless both sys- 
tems are described by the same low energy field theory 
([5]). Terms that break just the particle hole symmetry, 
such as d x <t>(d x 9) 2 , arc irrelevant at the HI-MI critical 
point, and do not change either the parity (cos (</>+)) 
or string (sin (</>+) ) order parameters in the insulating 
phases. What are the minimal perturbations that elimi- 
nate the sharp distinction between these two phases and 
gap out the critical point separating them? 

Consider the following perturbation of the Hamiltonian 
0: 



XP = H+ + X 



dx 



sin(2</>+) 



LL 



(ira) 2 
dx 

9 I ^jCOs(2</) + +x). (18) 



Here H LL is the critical Luttinger liquid theor y at the 
transition and x — arctan(A/<7i) and g = \fg\ + A 2 . 
Clearly P is a relevant perturbation at the critical point, 
leading to a gapped state even when the parameter g\ 
crosses zero, where the critical point would have been. 

Note that the operator P breaks both the particle hole 
and lattice inversion symmetries. This is seen for ex- 
ample in the representation of the single chain as two 
coupled fermionic chains (as in Ref. [l2|]). Written in 
terms of the fcrmions the perturbation is 



sin (2cf) + ) 



L.2 



H.c. 



(19) 



where ip R a and ip L for a = 1 , 2 are the right and left 
moving fcrmions in the two fermionic chains, respectively. 
(Note that the two chains labelled by a are the two fic- 
titious spin- 1/2 chains discussed in the paragraph pre- 
ceding Eq. ([5]).) This term is clearly odd under inver- 
sion (which changes ip R a < — > ^La) ana - particle-hole 
(tp R / L a < — > ^r/l a) transformations, and even under 
time-reversal. 

The perturbation P can be viewed as a local polariz- 
ing field that induces a small "dipole" moment on every 
lattice site. (By "dipole" we mean J xn(x)dx, i.e. the 
electric dipole moment that we would get if the particles 
where charged. It is not related to the real moment of 



the dipolar atoms.) It is clear from the second line of 
(fT8| that the term —g cos(20 + + x) acts to lock the field 



to the values 



-X/2 or cj) + w 7r — x/2. The den- 



sity concentration is shifted accordingly by ~ (x/27r)a 
away from the lattice points, resulting in a local dipole 
moment at these points. The shift of —x/2 also implies 
that both parity (Op) and string (0s) order parameters 
gain a finite expectation value in the presence of this per- 
turbation. Hence the sharp distinction between the Mott 
and Haldane insulating phases is lost. It is interesting 
to note in this regard the similar effect of a symmetry 
breaking field on a conventional quantum phase transi- 
tions involving spontaneous breaking of that symmetry. 
The field gaps out the critical point while inducing a 
non vanishing order parameter in the disordered phase. 
It is remarkable that breaking inversion symmetry has 
this effect in the HI-MI transition despite the fact that 
neither phase has broken lattice inversion symmetry (at 
least in an infinite system or in a finite system with pe- 
riodic boundary conditions). 35 

We note that the converse is also true: a finite ex- 
pectation value of both parity and string orders imme- 
diately entails broken inversion symmetry in the ground 
state. This is seen by writing the local symmetry break- 
ing field as a product of the two non local order param- 
eters sm(2(f) + ) = 2sin(0 + ) cos(0 + ). 

The connection between the string and parity correla- 
tions and the breaking of inversion symmetry can also be 
understood from the microscopic viewpoint sketched in 
Fig. [2] A non-zero expectation value of the parity oper- 
ator (Op) is associated with pairing of particle and hole 
fluctuations in the ground state, while (0s) ^ corre- 
sponds to alternate ordering of the particle and hole fluc- 
tuations. Having both implies organization of the parti- 
cle hole pairs in the form of ordered dipoles. A cartoon 
wave-function capturing this mixed phase is given by: 



=n[(i+d)bt+(i-d)b\ +1 



o: 



(20) 



Changing d between and 1 facilitates a continuous con- 
nection between the MI and HI phases. 

In an optical lattice system it is relatively easy to apply 
a perturbation that breaks lattice inversion symmetry. A 
second laser with double the wavelength can be used to 
produce a lattice of asymmetric double well o 33 ' 34 . Terms 
such as sin(20 + ) are then allowed by symmetry and will 
therefore be imminently generated in the effective field 
theory. 

An intriguing question that naturally arises is whether 
a phase described by a wave-function such as Eq. (|2"D)) 
can occur by spontaneous breaking of inversion symme- 
try. In fact, there is a natural route by which such sym- 
metry breaking can take place. As discussed in Sec. Ill A| 
the critical line between the MI and HI phases is stable 
when K + > \. For K + < \ sub- leading terms, such as 
cos(40 + ), become relevant and open a gap that would 
stabilize a new phase. In particular, if the coefficient 
of the cos(4</)_|_) term is positive, the 4>+ field is locked 
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to ±7r/4, which implies spontaneously broken inversion 
symmetry. In this case, the particles (or density max- 
ima) are effectively pinned to a point which is either 1/4 
lattice spacing to the left or to the right of lattice sites. 
What microscopic interactions are required to stabilize 
such a phase is an open question. It can only occur if 
the value of the Luttinger parameter K + goes below ^ 
before the transition to the DW terminates the critical 
line separating the MI and HI phases. This point will be 
investigated in a later study. 

A completely different type of perturbation that is dif- 
ficult to realize in cold atom systems, but is nevertheless 
worth mentioning, is one that breaks the U(l) symmetry 
associated with number conservation. Such a perturba- 
tion can be added to the continuum theory ([7|) as a term 
of the form A J dx cos(8 + ). In the spin chain case [Eq. 
0], this term corresponds to a transverse magnetic field. 
Clearly, this perturbation is relevant and would gap out 
the critical point. Inside the HI and MI this term pro- 
duces a finite density of kinks in the field + and therefore 
leads to vanishing of both expectation values (cos(0+)) 
and (sin(0 + )). Again the sharp distinction between the 
HI and MI phases is lost. 




U/2 (U+3Kt)/2 



FIG. 3: Evolution of the phase diagram with inter-chain cou- 
pling predicted from bosonization. (a) Plane of the nearest 
neighbor interaction V versus t± (inter-chain tunneling) at 
fixed on-site interaction U and no inter-chain repulsion, (b) 
In the plane of V versus Vj_ (inter-chain repulsion) at t± — 
and fixed U. 



D. Continuum limit of weakly coupled chains 

We move on to discuss a system of two coupled EBHM 
chains described by the microscopic Hamiltonian (|3j). 
The starting point for a low energy description at weak 
inter-chain coupling are the decoupled single chain field 
theories (|6]) which involve four bosonic fields, 4> a ,± (with 
a = 1,2 a chain index) and their canonical conjugate 
fields 9 a ,±. 

The continuum limit of the inter-chain Hamiltonian 
([2]) is now conveniently written in terms of symmet- 
ric and anti-symmetric combinations of the even fields 
of each chain: $± = (0i,+ ± 02,+) /v2, and 6± = 
± a ,+) /V§: 



, 96 ( / /7T^ \ „„„ /. , 



■ [cos, (v 7 ^- 



X COS (01,-) COS (02,- 



cos (y2§. 
•} 



J -97 

UX T~7T cos 



{■Kay v 
x [cos (0i,-) cos (9 2 ,-) + 



(21) 



(22) 



where the bare values of the parameters are g§ = Vj_a, 
.96 = Vj_aB- and gj = —tj_iraA-. (Here A, B are the non- 
universal constants that appear in Eq. (fTTj) .) We shall 
analyze how these terms affect the low energy physics of 
the decoupled chains © in the limit of weak t± and V± . 

The HI phase - The simplest case to consider is the 
case where the decoupled chains are in the HI phase. 



Because of the finite gap in each chain's spectrum, any 
inter-chain coupling term is irrelevant in this case, which 
implies that the chains remain essentially decoupled at 
weak coupling. Nevertheless, since the string order pa- 
rameter is a non local object in terms of either the density 
or the phase field it is destroyed by infinitesimal inter- 
chain tunneling. In appendix [Al we show that the string 
correlations decay exponentially with a correlation length 
that scales as (A/£j_) 2 , where A is the gap (see also Ref. 
[HI). It is easy to understand this effect by consider- 
ing the low energy theory ([6|) with the perturbation (|22|) . 
In the HI phase of a single chain the the field 0+ is lo- 
calized near ±7r/2, so that the string order parameter 
(sin(0 + )) is fixed around ±1. However, the perturbation 
t± cos (-\/26-), which hops a particle from one chain to 
the other, induces a kink in the string order parameter of 
each chain. Therefore at any finite coupling strength this 
perturbation destroys the string order by creating a finite 
density oc (t±/A) 2 of kinks in the Z2 order parameter. 

We note that the product of the string operators on the 
two chains remains invariant to the inter-chain tunneling 
because this perturbation always creates a kink on one 
chain and an anti-kink on the other. It is therefore tempt- 
ing to define the product Osl = sm (0+.i) sm (0+.2) as a 
generalized string order parameter for the ladder. Such 
an operator still has a non-zero expectation value in the 
HI phase when inter-chain hopping is introduced. How- 
ever, it gains a finite expectation value also in the MI 
phase in the presence of a finite tj_ and therefore cannot 
serve to distinguish these two phases. To see this we con- 
sider the decoupled chains in the MI phase 0+, Q ~ 0. The 
action of the inter-chain tunneling operator at point x 
produces a kink-anti kink pair corresponding to a change 
of 0+,i from to 7r at x and of 0+,2 from to — tt. At 
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the point x the values of the fields are <t>+,i ~ 7r/2 and 
4>+.\ ~ — tt/2, therefore such a point contributes —1 to 
the product string operator. We therefore expect that 
in the presence of finite coupling ij_ the product string 
operator would gain an expectation value even in the MI 
phase, which is of the order of the kink density (t±/ A) 2 . 

It is natural to ask what is the fate of the quantum 
phase transition from the MI to HI phase in the two chain 
system where we cannot define a clear cut distinction 
between the two phase. 

HI to MI transition - The critical theory for this transi- 
tion is a Luttinger liquid described by H + , Eq. © , with 
gi = 0. In this region, £T_ is gapped, with the cos (26 a _) 
term more relevant than the cos(2<^ a ,_) term. We may 
therefore assume that the # Qi _ fields are pinned at the 
minimum of the cosine potential. The effective low en- 
ergy Hamiltonian in the $± sector is written as follows: 



Hmi^hi = 



k± (d x e±f + 4- (d x $±f 



dx- 



9i 



(naY 



■ cos 



V2$_ 



cos 



+37 / dx 



(na)' 



■ cos 



v/20- 



where the bare values are u± = u± \ 1 ± 



(23) 



K+ = 



K + 1 J 1 ± 2 y^ + , and 177 oc t±. Here we have replaced 

cos(# Qi _) by its non-zero expectation value. Note also 

(2) 

that since the fields 6 a - are pinned, the g\,l term in Eq. 
(f2Tj) (which contains cos (</> ai _)) is strongly irrelevant at 
this critical point, and therefore it was omitted in Eq. 
(|23|) . We will, however, consider its effect on the phase 
diagram in what follows. 

The Hamiltonian (|23|) is identical to the effective 
Hamiltonian derived in Ref. [l4| for a system of two 
Bose-Hubbard chains coupled by an inter-chain hopping 
term. The only difference is that here, due to the ex- 
tended interaction terms, a wider regime of parameters is 
accessible in Eq. (|23|) . Of particular interest is the tran- 
sition point from the HI to the MI phase, where gi = 
(which was only possible in the non-interacting limit in 
Ref. [HI). 

For completeness, we will now review briefly the RG 
analysis of H MI _, HI in Eq. (S3]), following Ref Q- The 
leading order RG flow equations are: 



dgi 


H 


dl 




dg 7 




dl 


V 2K_ 


dK+ 


K%g\ 


dt 


16tt 2 


dk_ 


K\gl 4 


di 


16tt 2 




<J7 



8rr 2 



(24) 



The fate of the system is determined by the competition 
between the 37 and g\ terms, which contain the dual 
cosine terms cos (\/2@-) and cos (v / 2 < I > -), respectively. 
If <7i dominates, then the system will be in a Hl-likc phase 
(for g\ > 0) or an Mi-like phase (for g\ < 0). If g-j 
dominates, then a new phase is stabilized. Note that 
at the (71 — critical point, the 177 term is relevant for 

k_ > i. 

We can use the flow equations (j2"4")) to deter- 
mine the form of the phase boundaries that sepa- 
rate the new g t± dominated phase from the HI and 
MI phases. Imagine that we start from a point on 
the critical line between the MI or HI and the large 
g t± phases, (<7t_ L ,c,ffi,c)> where both g t±)C and g lfC 
are small. Integrating the leading order flow equa- 
tions, we find that this point is mapped to the point 

(t ±)C (Ao/A) 2 - 1 /( 2 *-), flliC (A /A) 2 -*+/ 2 -^-/ 2 ) ) where 

Ao/A > 1 is the RG scaling factor. Clearly, this point 
must also be on the critical line. Therefore the crit- 
ical line must be of the formi^ t± tC oc |5i, c | a where 

2-l/(2K + ) . 

a = - — - t~ — - ' Voy weak to intermediate mterac- 

2 — K + / 2 — 1\ _ / 2 

tions, this gives a > 1. 

To elucidate the nature of the new phase that forms be- 
tween the Mott and Haldanc insulators for non- vanishing 
t± , we consider the possible sub- leading interactions that 
may become relevant in this phase. The single chain even 
sector Hamiltonian H+ (Eq. JJJ) can contain higher har- 
monics of the form cos(2n^ Qj+ ) with n > 1. However, 
when expanded in terms of <i>± = (4>x,+ i 4>2,+) /v^j 
these terms are seen to be irrelevant since they always 
contain the odd field whose conjugate 0_ is pinned. 
The term cos (v / 8 ( E>+) can also be generated to higher 
orders in the bare couplings. (For example, such a term 
is generated during the RG flow to second order in g\.) 
This term has a scaling dimension of 2K + , therefore it is 
relevant for K + < 1. This represents a strong interaction, 
therefore one may still expect that for small to moder- 
ate U , V and inter-chain coupling, this term would be 
irrelevant. In that case, the intermediate finite t± phase 
is gapless, and characterized by power law correlations. 
Using the bosonized expression for the boson creation op- 
erator [Eq. (fTTj) ] , the off-diagonal correlation function in 
this phase is 



(b a (x)%(0)) 



/ e ie + (x)/V2 e - 



-(0)/V2\ 



f 



(25) 

where we have replaced cos (0^2-) and cos(-\/20-) by 
their non-zero expectation values. Therefore, K+ can be 
estimated numerically in the gapless phase by measuring 
the power law decay of this correlation function. We 
conclude that the phase penetrating between the HI and 
MI at finite t± is a SF phase which remains stable as long 
as the above correlation function decays with a power 
smaller than j-. 



In Ref. 1J] it was found that turning on inter chain 
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tunneling between two chains that are in the MI phase 
can drive a transition to a superfluid at a critical value of 
tj_ . Here we showed how this transition drops to t± = 
as one approaches the MI-HI critical point. 

The actual value of the Luttingcr parameter K + can be 
calculated reliably only in the weak coupling regime. To 
estimate it for strong coupling we use numerical DMRG 
simulations (Section llllj) . If either the intra-chain U, V 
or the inter-chain V± are increased sufficiently, K+ can 
be driven below the critical value of 1. Then the system 
would undergo a Kosterlitz-Thouless (KT)-type transi- 
tion to an insulating phase. 

We finally note that in the limit of vanishing inter- 
chain tunneling ij_ = 0, but finite interaction V± =/= 0, 
no intermediate SF phase is formed near the MI to HI 
phase boundary. The inter-chain coupling is in this case 
marginally irrelevant at the critical point, affecting only 
a renormalization of the Luttingcr parameter. 

HI to DW transition - This transition is of the Ising 
universality clas o 12 ' 24 . It is controlled by the odd sector 
single chain Hamiltonian H- [Eq. ©], which describes 
a competition between the terms cos(26L) and cos(20_). 
The transition to the DW phase occurs when the latter 
term becomes more relevant, which for small 172, 33 hap- 
pens at K_ = 1. For this value of K_ the theory ([5]) 
can be refermionized and written as a quadratic Hamil- 
tonian of Dirac fcrmions with a mass term proportional 
to gi and a pairing term proportional to 33 . The model is 
diagonalized when formulated in terms of two Majorana 
fields 2 ^. At the critical point 52 = 93 (and K_ = 1), one 
of these Majorana fields becomes massless. A masslcss 
Majorana theory is equivalent to an Ising critical point. 

H t± is clearly irrelevant at the critical point, since it 
contains cos (\/2S-) while the dual field <&_ is pinned. 
The only effect this term could have is a slight bending 
of the HI to DW phase boundary. 

On the other hand we show below that Hv ± is a rele- 
vant perturbation at the Ising critical point. First we can 
replace the operator cos(v / 2$_) — cos(v / 2$+) appearing 
in Hv ± by its non-zero expectation value. (Note that 
these operators involve only the even modes <p a ,+ of each 
chain, and the even sector Hamiltonian is not close to 
its critical point.) Now in order to determine the scal- 
ing dimension of Hy ± , we need the scaling dimension of 
the operators cos (0i/2,_) at the critical point. As men- 
tioned above, there is a correspondence, via the map- 
ping to a the Majorana theory, between these operators 
and the two Ising models that describe the odd sector 
Hamiltonain close to the critical point. The cos (4>i/2,-) 
operators correspond to products of the spin operators 
of the two Ising model a 23 ' 27 . However since only one of 
these Ising models becomes critical at the transition from 
HI to DW the operator cos (4>i/2,-) has the scaling di- 
mension of Ising spins in 1 + 1 dimensions at criticality, 
which is 1/8. Hy ± is thus strongly relevant. This can also 
be understood simply by the fact that at the transition 
to the DW phase, the density wave susceptibility of each 
chain diverges, and therefore the density-density coupling 



between the two chains is strongly relevant. Hy ± then 
drives a transition to a density ordered phase with the 
order parameters of the two chains locked to each other. 
The high V± phase is thus not distinct from the V± = 
DW phase. 

The HI-DW phase boundary at finite V± can be ob- 
tained using the same scaling argument as we used be- 
fore for the phase boundary between the HI or MI and the 
large t± phases. This gives that near the V± — critical 
point, the critical line is of the form Vj_ )C oc (V c — V) 7 ^ 4 , 
where V c is the critical value for the transition of the 
single chain from the HI to the DW phase (which is 
V c = (U + i-Kt) /2 in the weak coupling limit). 

In this section we studied weak coupling between two 
chains. The main results are summarized in the phase 
diagram plotted in Fig. [3l The case of an infinite array 
of chains is not expected to yield qualitatively different 
results in the weak inter-chain coupling limit. The main 
difference would be that the intermediate superfluid state 
would in this case support true long range order instead 
of power-law correlations. 



III. DMRG RESULTS 

In this section we present numerical results to support 
the predictions of the field theoretical analysis from the 
previous sections. To this end we compute the ground 
state and lowest excitations of the extended Bose Hub- 
bard model ([T]) and using the Density Matrix Renor- 
malization Group (DMRG)£L. The section is organized 
as follows. We first treat the case of a single unperturbed 
chain (Sec. IIII A[) at filling of n = 1 boson per site. Ex- 
panding on the results of Ref. Q , we focus on the inter- 
play of the non local string and parity correlations near 
the transition between the MI and HI phases. In Sec. 
IIIIBI we investigate how a perturbation that breaks the 
lattice inversion symmetry affects the MI-HI phase tran- 
sition. We find that such a perturbation gaps out the 
critical point and thus eliminates the phase transition, in 
agreement with the field theoretical prediction. In Sec. 
IIII CI we move on to address a two-leg ladder with only 
interaction coupling between the two chains. Finally in 
Sec. IIII Dl we analyze the two leg ladder with tunnel cou- 
pling. As expected we find that this coupling leads to 
an intermediate superfluid phase between the Mott and 
Haldanc insulators. 

Before proceeding let us give the essential technical de- 
tails of the numerical calculations. To expand the domain 
of the HI phase we add to (fT|) the next-nearest neighbor 
interaction (V/8) J2i Srii6ni+2- The DMRG calculations 
are performed with open boundary conditions, while 
keeping m = 250 states per block. As usual in bosonic 
problems, we also need to truncate the Fock space of site 
occupations. For calculations presented in this section 
we allow the four occupation states n = 0, 1, 2, 3. Includ- 
ing one more occupation state per site had a negligible 
effect on the results in a sample of representative calcula- 
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FIG. 4: (Color online.) (a) Phase diagram of a single unperturbed chain in the space of (U/t, V/t), obtained from the numerical 
calculations, (b) The "charge" (o) and "neutral" (□) gaps along the same cut through the three insulating phases. The charge 
gap vanishes at the MI-HI transition, whereas only the "neutral" gap vanishes at the transition to the DW phase, (c) The 
parity (o), string (□) and Density wave (o) order parameters as a function of V/t along a line of constant U/t = 6. The string 
and parity orders are defined as the square roots of (|12|l and (|13[) respectively. 
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FIG. 5: (Color online.) Decay of the parity and string corre- 
lations, and of the single particle density matrix at the HI-MI 
transition. The power-law fits are consistent with the field 
theoretical predictions for the relations between the different 
decay exponents (see text). 

tions. The maximum number of sites in the calculations 
was 256, including chains of length L = 256 and two leg 
ladders of L = 128. 



A. Unperturbed single chain 

To map the phase diagram of a single chain (Fig. Hla)) 
we compute the parity, string, and Density wave (DW) 
correlations in the ground state. The long distance be- 
havior of these correlations as a function of V for a par- 
ticular value of the on-site interaction U is plotted in Fig. 
0Jc). As expected, the MI is characterized by non van- 
ishing Parity "order", the HI by string order, and the 



DW phase by non decaying density correlations. The fi- 
nite value of the string and parity correlations seen in 
the figure at the MI-HI critical point is due to the finite 
system size (L = 256). The field theoretical analysis pre- 
dicts power-law decay of these correlation functions with 
a non universal power a, which is directly related to the 
decay exponent v of the single particle density matrix via 
the relation a = \/{Av). Along the line of critical points 
separating the HI and MI phases the exponents are pre- 
dicted to run in the range 1/4 < v < 1 . These predictions 
are consistent with the power-law fits of the relevant cor- 
relation functions at the critical point, as presented in 
Fig. [5] (See explanation below on how the critical point 
is located in the calculations). 

Also shown in Fig. 01a) is the superfluid (SF) phase 
at low U, V , which is identified by measuring the decay 
exponent of the single particle density matrix v. The SF 
phase is stable when v < 1/4. 

In addition to ground state correlations we compute 
the gap to "charged" and "neutral" excitations. The 
neutral excitation gap A„ is obtained by targeting the 
lowest excitation in the sector with exactly n = 1 par- 
ticles per site. The "charge" gap is defined by A c = 
E o (+1) + E o (-1)-2Eq{0), where E (±l) are the ground 
state energies of the system with one more/less parti- 
cle. There is an interesting complication in extracting 
the bulk gap in the HI phase. For open boundary condi- 
tions this phase supports low energy edge excitations. We 
can identify these states by inspecting the density pro- 
file of the wave functions. The appearance of the edge 
states coincides with the transition to the HI phase and 
facilitates the most precise determination of the transi- 
tion point in a finite system. In most cases, however, we 
are interested in the bulk properties. To extract the bulk 
"charge" and "neutral" gaps, we lift the edge excitations 
to high energy by applying a sufficiently strong field at 
the edges: V e d g e{Sni — Stil)- The gaps are plotted in Fig. 
0Jb) for a cut of the phase diagram at constant U/t = 6. 
It is seen that both the "charge" and "neutral" gaps van- 
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FIG. 6: (Color online.) Charge gap as a function of V/t 
at constant U/t = 6 for different values of the inversion- 
symmetry breaking perturbation. A is defined in (|26|) . 



ish at the transition from the MI to the HI phase. On 
the other hand only the "neutral" gap vanishes at the 
transition from the HI to the DW phase. 



B. Breaking of lattice inversion symmetry 



Following the predictions of section III CI we add the 
term 



5H = \Y J {5n i b\b i+ i + H.c.) 



(26) 
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to the Hamiltonian ([TJ) . This is one of the simplest terms 
that break the lattice inversion symmetry. Based on the 
predictions of III CI wc expect that the quantum critical 
point separating the HI and MI phases would be elimi- 
nated in the presence of this term and an adiabatic con- 
nection between the HI and MI would be facilitated. This 
is indeed what is seen from the calculated charge gap in 
the presence of the perturbation. In Fig. [5] we plot the 
charge gap along a cut through the phase diagram with 
U/t = 6, as a function of V/t. We see that for a non van- 
ishing value of A the gap does not vanish and the phase 
transition is eliminated. 



FIG. 7: (Color online.) Effect of inter-chain interaction cou- 
pling on the phase diagram, (a) Calculated DW order param- 
eter as a function of V/t at constant U/t — 6 for different 
values of the interchain interaction V± . (b) Calculated charge 
gap on the same cut through the phase diagram and the same 
values of Vj_. (c) Evolution of the phase boundaries with V± 
inferred from the calculations. 



values of V. The change in both phase boundaries with 
increasing V± are plotted in Fig. [7J The transition points 
are inferred from the vanishing of the "charge" gap at the 
MI-HI transition and the emergence of DW order. 



C. Two leg ladder with inter-chain repulsive 
interaction 

We move on to treat two leg ladders. The simplest 
coupling between two chains is via the density-density 
interaction V± J2i n i,i n 2,i- The field theoretical analysis 
predicted that this interaction is marginal at the transi- 
tion between the MI and HI phases. In other words it 
leaves the transition in tact, affecting only a renormaliza- 
tion of the Luttinger parameter which controls the decay 
of the single particle density matrix at the critical point. 

At the transition from HI to DW the interaction cou- 
pling V± is expected to shift the critical point to lower 



D. Two leg ladder with Inter-chain tunnel coupling 

Finally we include the inter-chain tunneling 
— 1± f&A fiB.i + H.c.J . According to the field theo- 
retical analysis of the previous sections this is a relevant 
coupling at the MI-HI critical point. An intermediate 
phase is predicted to occur between these two insulators. 
In fact the intermediate phase is expected to be a 
superfluid if the parameter K + > 1. It would seem 
that rather strong interactions are required to violate 
this criterion. However we note that K+ may be 
highly rcnormalized and not easily estimated from the 
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FIG. 8: (Color online.) Effect of inter-chain tunnel coupling 
on the phase diagram, (a) Calculated charge gap as a func- 
tion of V/t at constant U/t — 6 for different values of the 
inter-chain tunneling t±. The gaps were obtained by extrap- 
olation to an infinite chain through finite size scaling analysis. 
The fact that the estimate for the gap are sometimes slightly 
negative is due to errors in the extrapolation. Note the gap- 
less phase that appears between the two insulating phases, 
(b) Power-law fit for the spatial decay of the single particle 
density matrix for the same cuts through the phase diagram. 
The fact that v < 1/4 in the gapless phase confirms that it 
is indeed a superfluid. The transition to the gapped phase 
seems to occur at the universal exponent v = 1/4 consistent 
with a Kosterlitz-Thouless transition, (c) Evolution of the 
phase boundaries with t± inferred from the calculations. 



microscopic parameters. Therefore only the numerical 
calculations presented in this section can confirm the 
nature of the intermediate phase. 

The results of the calculations arc summarized in Fig. 
[5] The charge gap is shown as a function of V/t in a 
cut through the phase diagram at constant U/t. For 
increasing t± we see a that a gapless phase opens up 
between the MI and HI phases. Up to numerical accuracy 
and finite size effects the results are consistent with the 
prediction that the gapless phase is established for any 
non vanishing value of t± . The exponent v with which we 



fit the decay of the single particle correlation function is 
shown as a function of V/t. The fact that v < 1/4 in the 
domain of vanishing gap is consistent with a superfluid 
phase which is destabilized in a KT transition on crossing 
to either the HI or MI phases. We note that we find 
similar result also for different values of U/t for which 
the HI phase can be realized. 



IV. SUMMARY AND DISCUSSION 

The purpose of this work was to investigate fundamen- 
tal problems concerning the nature of non local string 
order parameters, which were brought into focus by the 
possible realization in systems of ultracold atoms£. In 
particular we addressed the question of how static per- 
turbations, which couple to local physical operators can 
nonetheless influence quantum phases and phase transi- 
tions that involve the highly non-local order parameter. 

A key result of this analysis is the discovery of a sur- 
prising connection between the string order and breaking 
of the lattice inversion symmetry (Note that breaking of 
particle-hole symmetry is also needed - but this is any- 
way broken at the outset in the Bose system) . In the con- 
text of this work the two phases that are characterized 
by a non local order parameter are the Haldane insula- 
tor, which sustains string order, and the Mott insulator, 
which supports non zero expectation of the parity op- 
erator. We find that a perturbation which breaks the 
lattice inversion symmetry gaps out the critical point be- 
tween the two phases. Correspondingly the sharp distinc- 
tion between them is eliminated, in the sense that non 
vanishing string order is induced in the Mott phase and 
parity expectation value in the Haldane phase. This is 
similar to the effect of a symmetry breaking perturbation 
on a conventional phase transition involving spontaneous 
breaking of the same symmetry. But curiously, neither 
the Mott or Haldane phase involve spontaneous breaking 
of the lattice inversion symmetry. 

The seemingly mysterious connection between lattice 
inversion symmetry and the string operators is elucidated 
by the effective long wavelength description of the prob- 
lem. The symmetry breaking field sin(20 + ) has the sym- 
metries of a local dipolc field: it is odd under both lattice 
inversion and particle-hole transformations. Incidentally, 
this operator can also be decomposed as a product of the 
long wavelength expressions for the string sin(0 + ) and 
parity cos(0+) operators. Thus in each phase only one 
(non local) factor of the local "dipole" field gains an ex- 
pectation value and so inversion symmetry remains in- 
tact. But this also implies that in each of these phases 
the external perturbation becomes in effect a direct cou- 
pling to the non local order parameter of the other phase. 
For example in the MI phase where (Op) = (cos0+) > 
we have that Asin20 + « X(Op)Os, which is essentially 
an ordering field for the string order. 

It is relatively easy to apply such a symmetry break- 
ing field in an optical lattice realization of the transi- 
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tion from Mott to Haldanc insulator. The principle has 
already been demonstrated successfully in experiments 
that created a lattice of asymmetric double wells using 
a secondary laser with half the wavelength of the main 
lase r 33 i 34 . A straight forward extension would be to os- 
cillate the secondary laser in time. This can be used to 
measure the dynamical response of this system to the lo- 
cal dipolc field, an interesting theoretical problem, which 
we leave for future work. 

The second part of our analysis addressed the effect of 
coupling between parallel one dimensional chains, which 
is inherent to realization with an optical lattice. As a 
simple model for the coupled chains we analyzed the case 
of two coupled chains, using a bosonization approach in 
section III Dl and numerical simulations with DMRG in 
sections IIII CI and IIII Dl The natural couplings to con- 
sider are inter chain repulsion due to the dipolar moment 
of the atoms or molecules and the inter-chain tunnel- 
ing. The interaction coupling turns out to be essentially 
trivial and does not change the structure of the phase 
diagram. The inter-chain tunnel coupling, on the other 
hand, has a dramatic effect. Like the inversion symmetry 
breaking perturbation discussed above, this perturbation 
eliminates the distinction between the MI and HI phases. 
However, instead of gapping out the critical point, it ex- 
pands it into a phase. More precisely, a weak inter chain 
tunneling gives rise to an intermediate supcrfluid phase, 
whose domain grows quickly with the coupling strength. 
The transition from the superfluid to either the Haldane 
or Mott insulator is Kosterlitz-Thouless like and occurs 
at the universal value of the decay exponent of the single 
particle density matrix. 

The evolution of the phase diagram with increasing 
inter-chain coupling for the case of many coupled chains 
and more general interactions deserves further theoreti- 
cal study. This may provide a controlled route for inves- 
tigating the crossover from one to two dimensions. Of 
particular interest is the question of possible generaliza- 
tions of the string orders to the case of two (or at least 
quasi one) dimensional systems. 
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APPENDIX A: STRING CORRELATIONS AT 
WEAK INTER-CHAIN COUPLING 

In this appendix we derive the effect of weak inter- 
chain coupling terms on the long range string correlations 
in the HI phase. For simplicity we will study this within 



the effective spin-1 model that results from restriction to 
three occupation states n = 0, 1, 2 on each sit o 21 ' 22 : 

(ij) * 

• l »" s ; /'»'<' 

(ij) « 

-M£]T( 5 r(^ + Sj)S+ +H.c.) 

(ij) 

~- ( s i s ? s I s j + sts*s*s~) (Al) 

where £; = \J~2 — 1 for a system single average occupa- 
tion. Here we shall set £ = 0, essentially neglecting terms 
that break the particle hole symmetry. Possible effects of 
these terms, which deserve further study, will be briefly 
discussed at the end of this section. 

Next we follow Kennedy and Tasaki and define a non 
local unitary transformation of the spins on all lattice 
sites 

U= ' e i7rS ^ (A2) 

jk\j<k 

The Hamiltonian gains an unusual, but nonetheless local 
form, in terms of the transformed spin variables: 

H = -jJ2SjSj+i-SjeM™$j+™Sj+i)Sj+i 

j 

-^E^I+i + yE^) 2 - (A3) 

j j 

H has an explicit Z2 x Z2 symmetry, generated by ir ro- 
tations of the transformed spins around the main axes. 
Furthermore the non local string correlations of the orig- 
inal spins map to standard two point correlations of the 
transformed spins. In the Haldane phase the Z2 X Z2 
symmetry is broken, and the transformed spin operators 
gain a finite expectation value. 

We can therefore treat the Haldane phase of H with a 
mean field approximation 13 , in which the four degenerate 
broken symmetry ground states are given by the product 
wavefunctions: 

1*0*) =n( cose l ) l + Csine|r ? l), t ). (A4) 

i 

The labels £,77 = ±1 correspond to the signs of the two 
Z 2 order parameters (Sf), and (Sf). The elementary 
excitations in this phase arc domain walls separating any 
two of the degenerate ground states. Interestingly such 
domain walls are generated by acting on a ground state 



13 



with the original spin operators: 
Sf = S?exp(i7rJ2§j) 

j>i 

Sf = cxp(i7r^5|)5f exp(iTr^S'J) 

j<i 

Sf = exp(^]Ts;)sf 



set by the average density and its harmonics k n — Imipo'. 



>l 



(A5) 



3<i 



Thus Sf creates a kink in 77 ((Sf)), Sf creates a kink in 
C ((Sf)), and Sf creates kinks in both order parameters. 
As in any Ising system, proliferation of a finite density 
of domain walls leads to destruction of the order. This 
is the mechanism by which inter-chain coupling destroys 
the long range string correlations. 

Consider first the perturbativc effect of the inter-chain 
interaction Hy ± = V± S^S^- The correction to the 
ground state to first order in V±/A consists of a single 
domain wall in the (Sf) order parameter in each chain. 
If we consider higher order terms, the weight of a config- 
uration with n domain walls per chain scales as (V±/A) n 
and the density of domain walls of (Sf) in the ground 
state is oc V±/ A. Therefore we expect the S x string cor- 
relations to decay exponentially over a correlation length 
£ x ~ A/Vj_. On the other hand the S z string order is 
left unmodified by this perturbation. 

We can apply a similar argument to the inter-chain 
tunneling term H t± = -t± Eii^fSJ + SfSf). 



APPENDIX B: ALTERNATIVE BOSONIZATION 
SCHEME 

In section |TI] we derived the continuum field theory in 
a rather indirect way. In the first step the EBHM ([I]) was 
mapped to the effective spin-1 model [AT| We then fol- 



lllGJ] for such 



lowed the procedure developed in Rcfs. 
spin systems which involves: (i) splitting a spin-1 chain 
into two effective spin-1/2 chains; (ii) fermionizing the 
parallel spin-1/2 chains with a Jordan- Wigner transfor- 
mation; (iii) Bosonizing the fermions. 

There is a standard procedure to "bosonize" a bosonic 
Hamiltonia n 17 ! 32 , and it is tempting to ask why not to 
use this more direct approach to derive the field theory 
which describes the transitions between the MI, HI and 
DW phases. This turns out to be not so trivial. 

The standard scheme is based on expansion of the fluc- 
tuations of the discrete particle density to slow modes 
describing small deviations from the fundamental period 



p(x) -> (p - -d x <p(xU 



a i(k n Xi — 2ri(p(xi )) 



(Bl) 



Thus, slow variations around the fourier component k = 
nko of the density are given by the field cos (2ncj)(x)). A 
finite expectation value for this field implies a static den- 
sity wave, with the periodicity A„ = 2vr/fc n = l/(np ) in 
the ground state. Clearly the period 2a density generated 
at unity filling by a strong nearest-neighbor interaction 
V is not captured in this expansion (here po = 1/a, with 
a the lattice spacing). It is therefore not surprising that 
the Haldane insulator, which involves fluctuations at the 
same length scale and is in some sense a precursor of the 
DW phase, cannot be described in this approach either. 

We now propose a modified bosonization procedure 
that will enable us to derive the low energy field the- 
ory ©. The key idea is to split a bosonic chain into 
two auxiliary chains without changing the total number 
of bosons. In other words, a single chain with one boson 
per site maps to a ladder with one boson per rung, or one 
boson in every two sites in each of the auxiliary chains. 
Next each of the half filled chains is bosonized separately 
in the standard way 



1/2 



E 



Sp a (xi) 



x > e 
mez 
1 1 
2a 7r 



d x 4> a (x i ) ) Y e 



(B2) 



where a = 1,2 is the auxiliary chain index. Note that 
the correct wave- vector to describe the DW has emerged 
from the inverse density of each of the split chains. We 
anticipate that the density wave of the physical chain 
will appear as in-phase locking of density waves in the 
two auxiliary half filled chains. 

We now map the EBHM ^ to a closely related model 
on the auxiliary ladder by taking b{ — > (bu + b 2 i)/ \/2 and 
rii — ► Tin + n2i- Choosing a slightly different extension 
of the BHM to the auxiliary ladder system should not 
change the essential structure of the phase diagram. The 
naive continuum limit of this model can now be taken 
by using the identities (|B2[) . For the on-site interaction 
term we obtain: 



U 



Ua 



dx 



r 1 



-) 2 - 4- cos (20+) cos ( 2 </ ) +) + -4( COS ( 2 0+) + cos(20_) N )l (B3) 
a 2 a z J - 



2 ^ v ' ' ' 2 

i 

where <j>± = cj)i zb 02, and we have kept only the most relevant terms. Similarly the nearest neighbor interaction term 
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leads to 
V 



/" r 1 1 

^(ni,i+n2,i-l)(ni,i+n2,j-l) « Va / dx —{d<j) + ) 2 — ^ 



cos(2c/> + ) cos(20_ ) ^ ( cos(20+) +cos(20_ (B4) 



Finally the hopping term translates to 

2 



i a, 0=1 



{a(de+f + a(<96L) 2 ) (1 + cos(26L)) - - cos(2</>+) cos(20_)(l + 2cos(20_)) 



s ,„ , , 2 4 ,_ , , 4 



47T 



(<90_) 2 + - cos(20+) + - cos(20_) - - cos(20_) 



(B5) 



Here 9± = (0i ± 6*2 )/2. Note that we have essentially the 
same degrees of freedom here as in section III Al In the 
MI and HI phases the operator cos(20_) is relevant and 
we can safely replace it with its expectation value C_ 
which is of order 1. In this case the above expressions 
simplify and give precisely the field theory ([6|) with the 
parameters: 



K+ = 



we expect 



also in the HI phase. 




,9i = .92 = Tr'a I - - V - 2tC 

i — + V + 2i(l + 2C*_) 




irat 



(B6) 



Note that the bare value of the Luttinger parameter K_ 
of the antisymmetric fields is large (K- > 2tt), consistent 
with taking C_ of order 1. 

Before closing this appendix let us make a few remarks 
concerning the new scheme. First note that the bare 
value of the Luttinger parameter K + implies a reasonable 
estimate for the phase boundary between the superfluid 
and insulating phases. The system should be superfluid 
when the renormalized Luttinger parameter K + exceeds 
2. Estimating this criterion with the bare value given 
above we get the approximate criterion (U + 2V)/t < 
7r 2 /4. By contrast the spin-1 mapping predicted that 
there is no superfluid phase in the region U > and 
V > 0. Indeed we do not expect the spin-1 mapping, 
which involves truncation to 3 occupation states, to hold 
when the on-site repulsion is not sufficiently strong. Thus 
the new scheme improves on the the spin one mapping in 
that it gives reasonable predictions for the phase diagram 
already at the level of the naive continuum limit. 

Another notable difference between the field theory of 
section Til Al and the new scheme is a ir shift in the defini- 
tion of the field (f> + . Hence in the bosonic scheme 0+ ~ tt 
corresponds to the MI phase whereas <j> + m to the HI. 
This is consistent with the physical interpretation of the 
variable <p + in the new scheme. Consider first the DW 
phase, which corresponds to a phase locked density wave 
on the two auxiliary chains. For this we need (f>i=(f>2 ~ 0. 
Since <t>+ is not critical at the transition from DW to HI 



APPENDIX C: THE STRING AND PARITY 
ORDER PARAMETERS 

The HI and MI phases do not support long range order 
in any local order parameter. Instead, they are charac- 
terized by the non-local string and parity order parame- 
ters (Eqs. (|15|14|) . respectively). In order to obtain the 
bosonized expressions for these order parameters in the 
effective field theory, we need to find their continuum 
limit. Here we propose a continuum form of these oper- 
ators using general considerations based on the asymp- 
totic "particle-hole" symmetry of the model at low en- 
ergy. More microscopic derivations can be found in Rcfs. 

In the following argument, we will assume that it is 
legitimate to truncate the Hilbcrt space to the three low- 
est occupation numbers (Srij = 0,±1). Sn can then be 
represented by a pseudospin-1 degree of freedom, defined 



by 5? 



The model (TTJ is then replaced by the effec- 



tive spin-1 model of Eq. ((4|, which leads to the same low 
energy effective theory (|7I8[) . We have neglected terms 
that break the symmetry between Sj = — 1 and 5*1 = 1 
( "particle- hole" symmetry). These terms appear in Eq. 
(|A1[) . They are irrelevant at low energies, so neglecting 
them should not change the long-distance behavior. 

Next, we need to find the continuum limit of the op- 
erators Op (j) and Os (j)- Taking the naive continuum 
limit of 6 P {j), we get 6 P - (since Y,j<i S j ~> 

j x dx' S z {x') = —cf)+(x), assuming that <fi + {— oo) = 
0) 28 ' 29 i 30 . To get a hermitian operator, e l< ^+ has to be 
symmetrized. To find the correct symmetrization, we 
note that Op (j) should be symmetric under a "particle- 
hole" transformation (which correspond under bosoniza- 
-> —</>+)• Therefore, Op should have the form 



tion to 



Op - A P cos {(j> + ) + 



(CI) 



where Ap is a non-universal constant, and we have 
truncated additional sub-leading operators of this sum. 
Os (j) also contains the same factor, but it should 
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be anti-symmetric under a particle-hole transformation, 
which suggests the general form 

6 S = A s sin (</)+) + B s d x cj) + cos(<£+) + ... (C2) 

The second contribution in Eq. (|C2[) is expected from 
naively bosonizing Eq. (fT5|) and taking Eq. (|C1[) into ac- 
count. However, a more careful treatment of the operator 
product expansion for this expression [26l(29l | shows that 
the first term is also present. The sin(0 + ) term is the 



most relevant one. From Eqs. (|C1IC2[) we see that in the 
MI phase, where </>+ is pinned around 0, we expect that 
(Op) 7^ 0, (Os) = 0, while in the HI + is pinned around 
±7r/2, therefore (6 P ) = 0, (6 S ) ^ 0. In order to test the 
validity of these results for the EBHM [Eq. |(TJ)], where 
particle- hole breaking terms exist, we evaluate the string 
and parity correlation functions in (| 1 211 3[) numerically 
across the MI— >HI transition. The results, summarized 
in Figs. I4|5( arc consistent with Eqs. (|C1IC2D . 
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